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Abstract 

We investigate an inequality constraining the energy and potential enstrophy flux spectra in two-layer and multi- 
layer quasi-geostrophic models. Its physical significance is that it can diagnose whether any given multi-layer model 
that allows co-existing downscale cascades of energy and potential enstrophy can allow the downscale energy flux to 
become large enough to yield a mixed energy spectrum where the dominant kr^ scaling is overtaken by a subdominant 
k~^^^ contribution beyond a transition wavenumber ki situated in the inertial range. The validity of the flux inequality 
implies that this scaling transition cannot occur within the inertial range, whereas a violation of the flux inequality 
beyond some wavenumber k, implies the existence of a scaling transition near that wavenumber. This flux inequality 
holds unconditionally in two-dimensional Navier-Stokes turbulence, however, it is far from obvious that it continues to 
hold in multi-layer quasi-geostrophic models, because the dissipation rate spectra for energy and potential enstrophy 
no longer relate in a trivial way, as in two-dimensional Navier-Stokes. We derive the general form of the energy and 
potential enstrophy dissipation rate spectra for a generalized symmetrically coupled multi-layer model. From this 
result, we prove that in a symmetrically coupled multi-layer quasi-geostrophic model, where the dissipation terms for 
each layer consist of the same Fourier-diagonal linear operator applied on the streamfunction field of only the same 
layer, the flux inequality continues to hold. It follows that a necessary condition to violate the flux inequality is the use 
of asymmetric dissipation where different operators are used on different layers. We explore dissipation asymmetry 
further in the context of a two-layer quasi-geostrophic model and derive upper bounds on the asymmetry that will 
allow the flux inequality to continue to hold. Asymmetry is introduced both via an extrapolated Ekman term, based on 
a 1980 model by Salmon, and via differential small-scale dissipation. The results given are mathematically rigorous 
and require no phenomenological assumptions about the inertial range. Sufficient conditions for violating the flux 
inequality, on the other hand, require phenomenological hypotheses, and will be explored in future work. 

Keywords: two-dimensional turbulence, quasi-geostrophic turbulence, two-layer quasi-geostrophic model, flux 
inequality 



1. Introduction 

It is now well-known that in two-dimensional Navier-Stokes turbulence, most of the energy tends to go towards 
larges scales and most of the enstrophy tends to go towards small scales, sometimes forming an upscale inverse energy 
cascade with energy spectrum scaling as k^^^^ and a downscale enstrophy cascade with k^^ scaling where k 

is the wavenumber Kraichnan [1] argued, diff'erently from Fj0rt0ft [4], that the direction of the two cascades can 
be justified via a thermodynamic argument in which we introduce, without proof, the assumption that the energy and 
enstrophy fluxes should tend to revert the energy spectrum from a cascade configuration to the absolute equilibrium 
configuration. The existence of forcing and dissipation arrests this tendency, thus keeping the system locked in a 
steady-state forced-dissipative configuration away from absolute equilibrium. 

Less well-known is the fact that there is a serious error with the original Fj0rt0ft argument: Fj0rt0ft claimed 
that the twin detailed conservation laws of energy and enstrophy alone imply that in every triad interaction group, 
more energy is transferred upscale than downscale. However, a more rigorous analysis shows that there exist triad 
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interaction groups in which more energy is sent downscale than upscale, and it is not obvious, without additional 
considerations, which group is dominant [5, 6]. Aside from this matter, the fundamental problem that underlies every 
other proof that utilizes only the twin conservation laws of enstrophy and energy, is that an additional assumption 
needs to be introduced to overcome the symmetry of the Euler equations under time reversal. Typical assumptions, 
such as the tendency of the energy spectrum to revert to absolute equilibrium, or the tendency of an energy peak to 
spread, typify ad hoc constraints imposed implicitly on the initial conditions that are needed to break the time reversal 
symmetry JtO- In Ref. 101 we counterproposed a very simple and mathematically rigorous proof that avoids the need 
for any ad hoc assumptions by considering the combined effect of the Navier-Stokes nonlinearity and the dissipation 
terms. The only assumption used by this proof is that the forcing spectrum is restricted to a finite interval [^i, ^2] of 
wavenumbers, however even that assumption can be relaxed to some extent, although not entirely eliminated ||8i Igfl . 

The essence of the argument in Ref. ff] is to show that for every wavenumber k not in the forcing range, the 
energy flux YlEik) and the enstrophy flux Yldk) satisfy the inequahty k^TlE(k)-ncik) < 0. Here, TlE{k) represents the 
amount of energy per unit volume transferred from the wavenumbers in the (0, k) interval to the wavenumbers in the 
(k, +00) interval, and llcik) is defined similarly for the enstrophy. From this inequality we then derive the following 
integral constraints for TlEik) and Yldk): 



These constraints imply a predominantly upscale transfer of energy and a predominantly downscale transfer of en- 
strophy. The original flux inequality k^TlEik) - Tlc{k) < itself can also be directly interpreted as a tight constraint 
on the downscale energy flux. 

The flux inequality is directly relevant to the cascade superposition hypothesis that was initially proposed in the 
context of two-dimensional Navier-Stokes turbulence \\V], according to which, for the case of finite small-scale 
dissipation viscosity, the downscale enstrophy cascade is accompanied with a hidden downscale energy cascade, 
associated with an accompanying small downscale energy flux. We stress that the existence of this small downscale 
energy flux is not in doubt. The hypothesis lies in the notion that it is part of a subdominant downscale energy cascade 
with both cascades contributing a k^^ and a k^^^^ term to the energy spectrum E(k), that are combined linearly, with 
similar linear combinations of terms to the generalized structure functions for all orders. This linear superposition 
principle can be proved for third-order velocity structure functions fl^]. A general consequence of this hypothesis 
is that, if the downscale energy flux associated with the k^^^^ term is strong enough, then a scaling transition in the 
energy spectrum from k^^ to k^^^^ should occur near a transition wavenumber k, x ^jr|u^ |euv, with t/j,,, the downscale 
enstrophy flux and e^, the downscale energy flux. The validity of the flux inequality for all wavenumbers k in the 
downscale inertial range of two-dimensional Navier-Stokes turbulence implies that the downscale energy flux is 
too weak to cause an observable scaling transition anywhere within the inertial range. On the other hand, it is far 
from obvious that the flux inequality will remain unconditionally valid in quasi-geostrophic models. A violation of 
the flux inequality beyond some wavenumber kt in quasi-geostrophic models would imply the occurrence of a scaling 
transition near that wavenumber A recent numerical simulation of a two-layer quasi-geostrophic model has indicated 
that a scaling transition is indeed possible [13]. Coexisting downscale cascades of energy and potential enstrophy have 
also been observed in more realistic models of atmospheric turbulence, the most recent being a primitive equations 
numerical simulation il4ll . 

The goal of the present paper is to extend the flux inequality to quasi-geostrophic models. We will specifically 
focus on vertical discretizations of the quasi-geostrophic model, namely the n-layer model, and the special case of the 
two-layer model, with all layers having the same thickness, in terms of pressure coordinates, on both models. From 
a physical standpoint, both models sacrifice the surface quasi-geostrophic dynamics at the bottom boundary, but they 
are otherwise good models of atmospheric turbulence for scales down to an estimated length scale of 100km 
I should like to emphasize from the beginning that in spite of any mathematical or phenomenological similarities, 
extending the flux inequality to quasi-geostrophic models is neither obvious nor straightforward. An overlooked 
fundamental difference between two-dimensional Navier-Stokes turbulence and quasi-geostrophic turbulence is that 
there are many more possible configurations for the dissipation terms in quasi-geostrophic models than there are in 
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two-dimensional Navier-Stokes. Dissipation terms are usually ignored because physical intuition alone may suggest 
that they should not have an effect on the nonlinear dynamics in inertial ranges. This line of reasoning ignores that 
the actual configuration of the dissipation terms can still have unexpected effects on the magnitude of the energy and 
potential enstrophy fluxes in the inertial range. These flux effects are the underlying matter of interest motivating the 
investigation initiated by the present paper 

The original motivation underlying the aforementioned numerical investigation lll3ll of the two-layer quasi-geostrophic 
model was to show that it can reproduce the Nastrom-Gage energy spectrum of the atmosphere lll6l4l9ll . However, the 
Nastrom-Gage controversy, reviewed to some extent in previous papers iflH^l, is not the main concern or motivation 
of this paper. Our main interest in this problem stems from the following considerations: first, quasi-geostrophic 
models are simple enough that they could be accessible to investigation via theoretical techniques developed for two- 



dimensional turbulence 1121 1 - 12611 . Furthermore, the possibility of being able to study a downscale energy cascade 
arising in the context of a two-dimensional model is particularly exciting from the point of view of the turbulence 
theorist, because it ties into the open question of why the downscale energy cascade of three-dimensional turbulence 
has intermittency corrections but the inverse energy cascade of two-dimensional turbulence does not f27, 28 1. Is it 
an the effect of dimension number or cascade direction? In light of such questions, an observable downscale energy 
cascade in a two-dimensional system is interesting in and of itself. 

Mathematical results concerning the flux inequality in quasi-geostrophic models can be organized into two cat- 
egories: (a) sufficient conditions for the satisfaction of the flux inequality within the entire inertial range; and (b) 
sufficient conditions for violating the flux inequality beyond some transition wavenumber kt within the inertial range. 
Results of the first type can be proved rigorously without ad hoc phenomenological assumptions on the behavior of 
the energy and potential enstrophy spectra. Results of the second type require the introduction of phenomenological 
assumptions about the distribution of energy and potential enstrophy between layers. Consequently, the scope of this 
paper has been limited to what we can prove rigorously. More powerful results that can be obtained by introducing 
phenomenological hypotheses will be explored in future publications. Because the details of our argument are very 
technical, we will now summarize the main argument of the paper as follows. 

For the generalized case of an n-layer model, we consider the general case of a streamfunction dissipation con- 
figuration, where for each layer the dissipation terms are given by a linear differential operator applied on the stream- 
function of the same layer, without entangling any streamfunctions of any other layers. The dissipation rate spectra 
for both energy and potential enstrophy are derived under this general configuration. Then, we specialize to the case 
of symmetric streamfunction dissipation, where we assume that the corresponding dissipation operators are identical 
layer-by-layer. We will show that under symmetric streamfunction dissipation the flux inequality is satisfied for all 
wavenumbers in the inertial and dissipation range. We note that this result is non-trivial since, beyond establishing 
cascade directions, it also implies bounds on the subdominant downscale energy flux, that are tight enough to keep 
the underlying downscale energy cascade hidden. For the case of the two-layer quasi-geostrophic model we consider 
an asymmetric configuration of dissipation terms and establish results of the form that if the asymmetry is sufficiently 
small, the flux inequality will remain valid. As was previously explained, we limit ourselves to results of this form 
because this is as far as one can go with rigorous proofs from first principles. 

From a physical standpoint, asymmetry in the dissipation between the two layers usually originates from the 
Ekman term, modeling the effect of friction with the surface boundary layer. However, for reasons that will be 
discussed more extensively at the conclusion of this paper, we will introduce an additional source of asymmetry via 
the small-scale dissipation terms by employing an increased viscosity or hyperviscosity coefficient at the bottom layer 
relative to the coefficient at the top layer. We believe that this asymmetric small-scale dissipation can facilitate a 
breakdown of the flux inequality, thereby allowing the downscale energy flux rate to be sufficiently strong to yield the 
transition to k'^^^ scaling in the inertial range. We will see that asymmetric small-scale dissipation indeed tightens the 
bounds on the parameter space wherein the flux inequality is satisfied. 

Another aspect of the dissipation term configuration, that will be shown to have significant impact on the flux 
inequality, concerns the modeling of the Ekman term. In a typical formulation of the two-layer quasi-geostrophic 
model, it is usually assumed that Ekman dissipation is dependent only on the streamfunction field of the bottom po- 
tential vorticity layer. However, an alternate formulation of the two-layer quasi-geostrophic model by Salmon ll29ll . 
requires that the Ekman term at the lower layer be dependent on the streamfunction fields of both layers. To explain 
why, one must recall that the two-layer model is an extreme vertical discretization of the full quasi-geostrophic model, 
which consists of a relative vorticity equation, a temperature equation, and additional constraining conditions. In a 
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general multi-layer model, the relative vorticity equations are discretized in horizontal layers that are interlaced with 
the discretization layers of the temperature equations. Thus, for the case of the two-layer model we have altogether 
5 physically relevant layers: the surface boundary layer corresponding approximately to lAtm, the lower relative 
vorticity layer at 0.75Atm, the temperature midlayer at O.SAtm, the upper relative vorticity layer at 0.25Atm, and the 
top boundary layer at OAtm. The potential vorticity equations are derived from the relative vorticity equations by 
eliminating the temperature field from the system of equations, thereby placing the potential vorticity field and the 



corresponding streamfunction field at the 0.25Atm and 0.75Atm layers. As noted by Ref. 02911 . the Ekman dissipa- 
tion term is dependent on the streamfunction field at the surface boundary layer near lAtm, which can be linearly 
extrapolated from the streamfunction field at the lower and upper layer (0.75Atm and 0.25Atm correspondingly). 
Consequently, even though the Ekman term is still placed on the lower-layer, owing to the linear extrapolation of the 
surface streamfunction field, it is dependent on the streamfunction field of both the lower and upper layers. 

It should be noted that for physical reasons, the potential vorticity layers need to remain fixed at 0.25Atm and 
0.75Atm respectively. This corresponds to the physical assumption that the two fluid layers have equal thickness, 
which is a necessary assumption for atmospheric modeling 1*3^. The surface layer driving Ekman dissipation, on the 
other hand, can be placed anywhere between the surface layer at lAtm and the lower streamfunction field layer at 
0.75Atm. When the surface layer and the lower streamfunction layer coincide, this corresponds to the usual standard 
Ekman term. When the two layers do not coincide, it corresponds to the more general case of extrapolated Ekman 
dissipation. For the present paper, we retain generality by parameterizing the placement of the surface boundary layer 
via an adjustable parameter//, and show that our main propositions are valid for the entire range of the parameter 
We will see that an increasing separation between the surface layer and the bottom potential vorticity layer tightens the 
bounds on the parameter space wherein the flux inequality is satisfied. For oceanographic modeling, as well as for the 
purpose of satisfying basic scientific curiosity, it would be interesting to consider two-layer quasi-geostrophic models 
with layers having unequal thickness. Due to mathematical complications, we will not pursue this generalization in 
the present paper Nevertheless, the importance of symmetric vs. asymmetric Ekman dissipation in the context of 
oceanographic modelling is a relevant problem that has been investigated by a previous study llsill . 

Admittedly, both Salmon's idea of extrapolated Ekman dissipation and my idea of diff'erential small-scale dis- 
sipation can be considered controversial. On the other hand, in the context of investigating the flux inequality, it 
is important to be thorough about considering every interesting configuration of the dissipation terms, to determine 
how much impact various choices of dissipation term configurations have on the robustness of the flux inequality. 
Furthermore, as will become apparent from the results of this paper, the dissipation configurations explored here are 
good candidates for a dissipation filter that could violate the flux inequality and ensure a controlled downscale energy 
dissipation rate in numerical simulations that exceeds the restrictions that are typical in two-dimensional turbulence. 

The paper is organized as follows. In section 2 we give the governing equations for the generalized multi-layer 
model and discuss its conservation laws, the definition of the energy spectrum E(k), potential enstrophy spectrum 
G(k), and their relationship via the streamfunction spectrum Cap{k). In section 3, after a brief recapitulation of the 
flux inequality for the simple case of two-dimensional Navier-Stokes turbulence, we establish the flux inequality 
for a generalized multi-layer quasi-geostrophic model under symmetric streamfunction dissipation. In section 4, we 
consider asymmetric dissipation configurations for the special case of a two-layer quasi-geostrophic model, where we 
derive various sufficient conditions for satisfying the flux inequality. Conclusions and a brief discussion are given in 
section 5. 

2. The generalized multilayer model and conservation laws 



Following my previous paper 112011 . we write the governing equations for the generalized multi-layer model in 
matrix form: 

dqa 

— + J{^a, qa) ^da+ fa, (3) 
Ot 

P 

Here i{/a represents the streamfunction at the a-layer, qa represents the potential vorticity at the a-layer. Slap is a linear 
operator encapsulating the dissipation terms, and /„ is the forcing term acting on the (T-layer The index a takes the 
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values a - 1 , 2, . . . , « representing the layer number, for a model involving n layers. Sums over indices, such as in the 
sum over the index fi in the dissipation terms above, are assumed to run over all layers 1 , 2, . . . , n, unless we indicate 
otherwise. It is also assumed that the streamfunction and the potential vorticity are related via a linear operator 
according to: 

qa{^, = ^ ^afi^fii^, t). (5) 

P 

The above equations encompass both the two-layer quasi-geostrophic model and the multilayer quasi-geostrophic 
model, on the assumption that we neglect the y6-effect, arising from the latitudinal dependence of the Coriolis pseud- 
oforce. This is a reasonable assumption for Earth, especially if we restrict our interest to a thin strip of the Earth's 
surface, oriented parallel to the equator Baroclinicity instability is accounted for by the forcing term /„, and implicit 
in the entire argument is the assumption that it forces the system at large scales only. This assumption, originally pro- 
posed by Salmon 1 2^ 32 ], is the only physical assumption implicit in the theoretical framework of the flux inequality. 



and it has been corroborated numerically Ill3[ 13311 . 



For the sake of simplifying our analysis, we assume that all fields are defined in an infinite two-dimensional 
domain. Then we can write the Fourier expansions for the streamfunction ij/ and the potential vorticity q as follows: 



if/ai^, t) - I lAffC^, t) exp(/k ■ x) dk, (6) 
^Q,(x, f) - I §Q.(k, f) exp(/k ■ x) dk. (7) 



We assume that the operator is diagonal in Fourier space. This means that the relation between the streamfunction 
and the potential vorticity, in Fourier space, reads: 

§„(k,f) = ^L.MIIk|l)'/'<.(k,0. (8) 

P 

Here ||k|| represents the 2-norm of the vector k. We also assume that is symmetric with - ■^pa- This 
implies that Lap{k) - LpaiK) for all wavenumbers k. For quasi-geostrophic models, the matrix La/sik) is non-singular 
for all wavenumbers k > Q, due to being diagonally dominant, and we assume that to be the case in our abstract 
formulation given above. Consequently, there is an inverse matrix L^j^ik) which defines the inverse operator . To 
accommodate a possible singularity at ^ = we assume that at wavenumber A; = 0, in Fourier space, the corresponding 
field component is for all fields. This is equivalent to subtracting the mean field and considering only the field 
fluctuation around the mean. 



2.1. Conservation laws 

We will now show that the generalized layer model, in the absence of dissipation, conserves the total energy E 
and the total potential enstrophy G under very general conditions on the operator For any arbitrary scalar field 
f(x, y) we write the corresponding volume integral using the following notation: 

«/»= rr f(x,y)dxdy. (9) 

J JR- 

We define the total energy E over all layers, and the layer-by-layer total potential enstrophy Ga for layer a, as E - 
- Yjaii'l^aQa)) and Ga - ((^^)). The purpose of the minus sign in our definition of the total energy E is to maintain 
consistency with the notation and sign conventions used by my previous paper 1.20.1 . Specifically, we will show that 
the potential enstrophy is conserved on a layer-by-layer basis unconditionally regardless of the details of the operator 
^ap- Conservation of the total energy E, over all layers, on the other hand, requires that the operator Jifap be symmetric 
and self-adjoint. By symmetric we mean that the operator satisfies = -^pa- To define the self-adjoint property, 
consider two arbitrary two-dimensional scalar fields f(x,y) and g{x,y). We require that every component of the 
operator must satisfy {{fi^apg))) = €i-^apf)8)) for any two fields f(x,y) and gix,y). This self-adjoint property, 
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so defined, follows as an immediate consequence of our previous assumption that the operator is diagonal in 
Fourier space. In the proof given below, however, there is no need to use the stronger assumption of diagonality. 

The proof is based on the following properties of the nonlinear Jacobian term. If a(x, y) and h{x, y) are two- 
dimensional scalar fields that satisfy a homogeneous (Dirichlet or Neumann) boundary condition, then we can show 
that ((7(fl, /?))) = 0, using integration by parts. Then, we note that, as an immediate consequence of the product rule of 
differentiation, given three two-dimensional scalar fields a{x,y), b{x,y), and c(x,y) we have 

{{J(ab, c)» = {{aJ(b, c)» + ibJia, c))) = 0, (10) 

from which we obtain the identity 

iaJ{b, c))) = {{bJ{c, a))) = {{cJ{a, b))). (11) 

Now, let us go ahead and drop the dissipation and forcing terms and write the time-derivative of the potential vorticity 
qa as qa - -JiiJ/a, qa)- Then, the time derivative of the streamfunction if/a reads: 

Differentiating the total potential enstrophy Ga for the a layer with respect to time and employing the identity 
given by Eq. ( fTTT i immediately gives: 

Ga = liqaqa)) = -2{{qaJ(l//a,qa)}) = -2«(/r„y(^<,, g'^-))) = 0. (13) 

Here, we note that from the definition of the Jacobian J(qa, qa) - 0. This establishes the layer-by-layer conservation 
law of potential enstrophy, unconditionally, as claimed. To show the energy conservation law, we differentiate the 
total energy E with respect to time and obtain: 

E = -(d/dO = - - Yjii^'r^^)) (14) 

a a a 

= J]fc^„^'y(iA/J, qp)}} + J^UaJiifTa, qa)}} (15) 

iY/3 a 

= qi^)^Kfiqa» + Yu^{qaJ{4>a, <AJ» (16) 

aft a 

= ^«/(^/3, qp)J^p^qa)) = Y^Ji'Pp, qp)'l^p}} (17) 

= ^«7((A/?,-A/?k/j» = 0. (18) 



Note that the self-adjoint property is applied at Eq. (fT&t . and the symmetric property is applied at Eq. (ITTt . This 
concludes the proof. 

2.2. Definition of spectra 



Following Frisch Il34ll . we define spectra for the energy and potential enstrophy using the bracket notation intro- 



duced in my previous paper f^O*]. Consider, in general, two arbitrary two-dimensional scalar fields a{x) and b{x). 
Let a^''{x) and b^'^ix) be the fields obtained from a(x) and bix) by setting to zero, in Fourier space, the components 
corresponding to wavenumbers whose norm is greater than k. Formally, a^''(x) is defined as 



Jr2 Jr2 



a^''(x)^ I dxo I dko exp(/ko ■ (x - xo))fl(xo), (19) 



47r2 

with H(x) the Heaviside function, defined as the integral of a delta function: 



if if X e (0, +oa) 

H{x)= \ 6{T)dT^\ 111, ififx = . (20) 

if if jc € (-00,0) 
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Obviously, b^''{x) is defined similarly. We now use the two filtered fields a^''(x) and b^''{x) to define the bracket (a, b)i^ 



as: 

{a. 



b), ^ j^ dx {a<''(x)b<\x)) (21) 
= - r dQ(A) {[a(kAe)b(kAe) + a(kAe)b*(kAe)]) . (22) 

2 jAeS0(2) 

Here, a(k) and £(k) are the Fourier transforms of a{x) and b(x), SO(2) is the set of all non-reflecting rotation matrices 
in two dimensions, dD.{A) is the measure of a spherical integral, e is a two-dimensional unit vector, and (•) represents 
taking an ensemble average. The star superscript represents taking the complex conjugate. Note that Eq. (ISTT i is the 
definition of the bracket, and Eq. (l22l i follows from Eq. (|211 as a consequence. The bracket satisfies the following 
properties: 

{a,b}k = {b,a}i,, (23) 
(a, b + c)i^ - (a, b)i^ + {a, c)i; , (24) 
(fl + b, c)k = (a, c)i + {b, c)^ . (25) 

Moreover, every (Q'y6)-component of the operator ^a/i is self-adjoint with respect to the bracket, which gives 

(^^a/^a, b^^ = (fl, .^apb)^ = Lai}(k) {a, b)i, , (26) 

and the same property is also satisfied by every component of the inverse operator : 

{■^a'pa, b)^ = {a, ^^^b)^ = L-lik) {a, b\ . (27) 

Using the bracket, we define the energy spectrum E(k) - - J]„ {il^a, qa}k^ ^nd we also define the layer-by-layer 
potential enstrophy spectrum Gaik) = {qa, qa)k ™d the total potential enstrophy spectrum G(k) - 2„ Ga{k). Unlike 
the case of two-dimensional Navier-Stokes, where the enstrophy and energy spectra G(k) and E(k) are related via 
a simple equation, G{k) - k^E{k), in the generaUzed layer model, the potential enstrophy spectrum and the energy 
spectrum are related indirectly, as shown below: 

Define the streamfunction spectrum Capik) - (i//a, Then, via the properties of the bracket above, the energy 
spectrum E(k) reads 



E(k) = - ^ q,}, = - ^ Ua, 2 ^afiM = - 2 L„fi(k) -A/i), (28) 
^-Y,Laii(k)Cap(k), (29) 



and the potential enstrophy spectrum Ga(k) reads 



G(k) = i^'^k = Z (Z ^"f'^"^ Z ^"y'^yj ^^^^ 

= Lafiik) Y ^ayl^y) = ^ ^«/?(^)'L-r(^) ('A/*^ "Ay), (31) 

ffjS ^7 ' h ally 



Yi^al3(k)Lay(k)C^y(k). (32) 



al3y 



Thus, they are related only indirectly via the streamfunction spectrum Cap{k). 

We note that for a p, Capik) may take positive or negative values. For the case a = /3 we define Ua{k) = 
{ilfa,^/^a)k^ which is always positive (i.e., Ua{k) > 0). Then we note that 2\Ca0(k)\ < Ua{k) + Up(k). We can use this 
inequality to show that if the matrix Lapik) satisfies the diagonal dominance condition 



Upik) > 0, for ai^p. 
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(33) 



(34) 



then the energy spectram E{k) is always positive. We give the proof in Appendix A Both the two-layer quasi- 
geostrophic model and the multi-layer quasi-geostrophic model satisfy this diagonal dominance condition. As for 
the layer-by-layer potential enstrophy spectra Ga{k), it is immediately obvious that they are unconditionally always 
positive, regardless of the form of the matrix Lapik), since by definition Ga{k) - (qa, qa)k- 



3. Flux inequality for the n-layer model 

We now turn to the main issue of identifying sufficient conditions for satisfying the flux inequality k^YlEik) - 
^cik) < for quasi-geostrophic models. Let us recall that the energy flux spectrum Y\-E{k) is defined as the amount 
of energy transferred from the (0, k) interval to the {k, +oo) interval per unit time and per unit volume. Likewise, the 
potential enstrophy flux spectrum Ylcik) is the amount of potential enstrophy transferred from the (0, k) interval to the 
{k, +oo) interval, again per unit time and volume. Assuming a forced-dissipative configuration at steady state and that 
the wavenumber k is not in the forcing range, the energy and potential enstrophy transferred into the (k, +oo) interval 
eventually are dissipated somewhere in that interval. It follows that we may write the flux spectra I\-E{k) and licik) as 
integrals of the energy and potential enstrophy dissipation rate spectra DE(k) and Ddk): 

DE(q)dq, (35) 

Dc(q)dq, (36) 

which implies that 

k^UE(k)-Uc(k)^ [k^DE{q)-DG(q)]dq^ A(k,q)dq. (37) 

Jk Jk 

We see that a sufficient condition for establishing the flux inequality is to show that A(k, q) < for all wavenumbers 
k < q. It is also easy to see that A(k, q) > for all wavenumbers k, < k < q is sufficient for establishing the violation 
of the flux inequality for all wavenumbers k > k,. 

For the case of two-dimensional Navier-Stokes turbulence, the dissipation rate spectra DE{k) and Dc{k) are related 
via Dc{k) - k^DEik). This immediately gives A(^, q) = k?DE{q) - Ddq) = (k^ - q^)DE{q) < for all wavenumbers 
k < q (since DE{k) > 0), which in turn gives the flux inequality k^'Il-Eik) - ^cik) < 0. The physical interpretation of 
this inequality is that when we stretch the separation of scales in the downscale range, the energy dissipation rate at 
small-scales vanishes rapidly. As a result, most of the injected energy cannot cascade downscale although, as noted 
previously ifloilTlll . a small amount of energy is able to do so. As we have seen in the previous section, for the case 
of quasi-geostrophic models, the energy and potential enstrophy dissipation rate spectra no longer have a direct and 
simple relation with each other, so the validity of the flux inequality needs to be carefully re-examined. 

For the general multi-layer quasi-geostrophic model, the relationship between the potential vorticities qa and the 
streamfunctions [j/a is given by 

qi = V^iAi +jJikl(if/2 - (Ai), 

qa = V^l/r^, - Aakliif/a - ^a-\) + l^ak^i^a+l - >Pa), for 1 < < «, 

q„ = V^i/r,, - /l„^(i/'„ - iA«-l)- 

Here, kji is ffie Rossby wavenumber and /!„ and are the non-dimensional Froude numbers, given by 

\hi p2-p\ ^ . 
Aa = -■; , for 1 < a < n, 

2 ha Pa -Pa-1 

\hi p2-p\ . . 
jia , for \ < a < n, 

2 haPa+l - Pa 
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with Pa the average density of layer a, and ha the average height of layer a (in pressure coordinates). The definition of 



the non-dimensional Froude numbers was adjusted with a 1 /2 numerical factor, from the one given by Evensen 113511 . 
to ensure agreement with the formulation of the two-layer quasi-geostrophic model given by Salmon |29] for the case 
n - 2. The components of the corresponding matrix Lap(k) are given by 

{-k^ - i-i\k^, if Of = 1 

-k^ - (Aa + iJ.a)kl, if 1 < a < « 
-k^ - A„kj^, \f a - n, 

La,a+\{k) - iiak%, for 1 < a < ti, 

La,a-i(k) — Aak\, for 1 < a <n. 

In the present paper we limit ourselves to the special case of a symmetrically coupled multi-layer quasi-geostrophic 
model, where we assume that the layer thickness ha is the same for all layers, thereby yielding a symmetric matrix 
Lafiik) such that La^a+i{k) = La+\^aik) for all 1 < a < n. 

To consider the flux inequality for this general «-layer model, we begin with writing the dissipation rates Dgik) 
and Dcik) for the energy and potential enstrophy in terms of the streamfunction spectrum Capik). We assume that the 
dissipation operation 2>ap is diagonal in Fourier space and that the Fourier expansion of the dissipation term Siap^'p 
reads: 

{%pil^p){x,t) ^ r £)„;,(||k||)iA/j(k,f)exp(/k-x)dk. (38) 



Then, in Appendix B we show that the energy dissipation rate spectrum Dsik) and the layer-by-layer potential en- 



strophy dissipation rate spectra Dc„ (k) are given by 

DE(k) = 2 2 Dap{k)Cap{k), (39) 

Dc„ {k) = -2 2 Lap{k)Day{.k)Cpy{k). (40) 

Note that in order for the dissipation terms to be truly dissipative, the dissipation spectra Dsik) and Dc(k) need to be 
both always positive for all wavenumbers k. From the general form of the above equations this is not readily obvious. 
However, for simpler configurations of the dissipation operators, the above expressions for DE{k) and Dc{k) simplify 
considerably, thereby making it possible to establish that they are both always positive. These expressions also 
underscore the main difference between two-dimensional Navier-Stokes turbulence and quasi-geostrophic turbulence 
and the reason why the flux inequality becomes a non-trivial problem in the later case. Unlike two-dimensional 
turbulence, and in spite of the twin conservation laws of energy and potential enstrophy, the dissipation rates Dgik) 
and Dcik) are no longer related by any simple relation of the form Dc{k) = k^Dgik). 

We restrict our attention to the case where the dissipation operators at every layer involve only the streamfunction 
of the corresponding layer, with no explicit interlayer terms. This can be arranged in terms of a linear operator 
applied to the streamfunction iffa- If Da{k) is the spectrum of the positive-definite operator then for the case of a 
dissipation term da - ^a^a, we have Dapik) - SapDp(k), with dap given by 

^"^ = {o, iia/p ■ ^41) 
We designate this case as streamfunction-dissipation. The DE{k) and Dc^{k) simplify as: 

DE(k) = 2 2 Dal3{k)Cap{k) = 2 ^ 6apDp{k)Cal3{k) = 2 ^ Da(k)Caaik) = 2 ^ Da(k)Ua(k), (42) 

ap ap a a 

Dc^ (k) = -2 ^ Laii{k)Day{k)Cpy(k) = -2 2 Lap{k)6ayDy{k)Cpy(k) = -2 ^ Laii{k)Dp{k)Cal}{k). (43) 



Py J3y 
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Note that for D„{k) > 0, it follows that DE{k) > 0, but it is not obvious that the same result extends to Dc^{k). 
However, if we further assume that the same operator is used for all layers, i.e. D„{k) = D{k), then we have the more 
specialized case of symmetric streamfunction-dissipation, and the dissipation rate spectra Dsik) and Ddk) can be 
simplified further to give: 

Dsik) = 2 2 D,{k)Ua{k) = 2D{k) ^ Ua{k) = 2D(k)U(k), (44) 

a a 

Dcik) = Yj ^G.(fc) = -2 2 Lai}(k)Di3(k)Cai3(k) = 2D(k) - ^ Up(k)Cap{k) = 2Dik)E{k). (45) 

a afi V afi . 

Now, D{k) > implies both Dsik) > and Ddk) > 0. 

It follows that, under symmetric streamfunction dissipation, A(fe, q) is given by 

Mk, q) = k'DEiq) - Dciq) = k^D(q)U(q) - D(q)E(q) = D(q)[k^U(q) - E(q)], (46) 

and since D{q) > 0, the validity of the flux inequality is dependent on the sign of the factor k^Uiq) - E(q). That sign 
is in turn intimately related with the expression jaik, q) defined as; 

ya{k,q)^k^ + Y^L„p{q). (47) 

p 

Note that for the case of two-dimensional Navier-Stokes, L{q) becomes a 1 x 1 matrix with L\ i iq) - q^, thus jaik, q) - 
k^ - q^, which is negative when k < q. For more generalized n-layer quasi-geostrophic models, the expression y„(A;, q) 
continues to be given by jaik, q) - k^ - q^ which remains negative when k < q for all layers a. We will now show 
that: 

Proposition 1. In a generalized n-layer model, under symmetric streamfunction dissipation da — +&4'a with spec- 
trum D(k) > 0, we assume that La/siq) > when a ^ f3, and Lapiq) — Lpaiq), and jaik, q) < when k < qfor all a. It 
follows that: 

A(k, q) < Diq) ^ Jaik, q)Ua(q) < 0. 

a 

Proof. We begin by recalling from [Appendix A that E{q) can be rewritten as 

E(q) = - ^ La^iqWaiq) - ^ ^ Lap(q)[2Cap(q) - Uaiq) - Up(q)] . (48) 

a/J a/3 

It follows that k^U(q) - E(q) satisfies: 

k^Uiq) - E(q) ^k'Yj Uaiq) + Yj La/siqWaiq) + ^ ^ Lapiq)[2Cap(q) - Uaiq) - Upiq)] (49) 

a ap lyP 

<k^Y Ua(q) + Yj LapiqWaiq) = + ^ L„^(^))f/„(^) (50) 

a cy/3 a 

a 

The inequality uses the assumption Lapiq) > combined with the triangle inequality 2Caj3(q) ^ Uaiq) + Upiq) of the 
streamfunction spectra. It follows that 

Mk, q) = Diq)[k^U(q) - Eiq)] < D(q) ^ ya(k, q)Ua(q) < 0, (52) 

a 

since D(q) > 0, Uaiq) ^ 0, and jaik, q) < 0, thereby concluding the proof. □ 
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The above result establishes the unconditional validity of the flux inequality for generalized «-layer quasi-geostrophic 
models under symmetric streamfunction dissipation. We note that the condition Lapiq) > is needed to establish that 
the energy spectrum E{k) is always positive, and all physically relevant quasi-geostrophic models will also satisfy the 
condition jaik, q) < for all k < q. As we have already argued, for any general n-layer quasi-geostrophic model, 
we have yaik, q) - k^ - (f- for all layers a, so the assumption is mathematical and does not impose any physical 
constraints in the model's formulation. No other restrictions are needed by the above proposition. Physically, this 
means that under symmetric streamfunction dissipation, the behavior of any generalized n-layer model will be similar 
to two-dimensional turbulence, where the subdominant downscale energy cascade is too weak to cause a transition 
from k'^ scaling to k^^^^ scaling in the downscale inertial range. 



4. Flux inequality in a two-layer model 

The previous results, derived over a general n-layer quasi-geostrophic model also apply to the special case of a two- 
layer quasi-geostrophic model. Consequently, the flux inequality will be satisfied by any two-layer quasi-geostrophic 
models under symmetric streamfunction dissipation. We will now concentrate on investigating the validity of the flux 
inequality in two-layer quasi-geostrophic models with asymmetric dissipation. 

4.1. Model formulation 

The two-layer quasi-geostrophic model can be formulated in terms of two potential vorticity equations of the form 
= (53) 

ot 

^+m2,q2)^fi + d2, (54) 
ot 

with the relationship between the potential vorticities qi, q2 and the streamfunctions i^i, t//2 given by 

?i = VVi +y(-A2-'Ai), (55) 

?2 = vV2-y(<A2-'Ai)- (56) 

Here qi, if/i correspond to the top layer and q2, 4^2 correspond to the bottom layer. As explained in the introduction, 
we situate the top layer at pi = 0.25Atm and the bottom layer at p2 = 0.75Atm. In terms of the generalized layer 
model, Eq. i55[ and Eq. ( |56] | correspond to an operator with spectrum Lapik) given by 



m 



a(k) b(k) 
b(k) a(k) 



(57) 



with a{k) and b(k) given by a(k) - k^ + k^fl and b{k) - -k^. Using diff^erential hyperdifFusion at the small scales and 
extrapolated Ekman dissipation at the bottom layer gives 

t/i = v(-l)'^+iv2"+Vi. (58) 
d2 = (v + Av)(-1)'^+'V2^'+V2 - VfiVV.. (59) 

Here we assume that the hyperdifFusion is stronger at the lower layer, thus Av > 0. Furthermore, the Ekman term is 
given in terms of the streamfunction i^j at the surface layer (ps = lAtm) which is hnearly extrapolated from ij/i and 
1/^2 and it is given by iffg - A^2 + jU/li/'i, with A and ji given by 

- and ^ . (60) 



P2-P\ Ps-P\ 

In other words, i/tj is defined so that, plotted on apressure-streamfunctionplane, the three points (p,, ijjs), {p\,il/\), (p2, 1^2) 
are geometrically collinear. Using p\ - 0.25Atm and p2 - 0.75Atm gives A - 3/2 and /i = -1/3. It is worth not- 
ing that for any arbitrary placement of the top and bottom layer that satisfies Q < p\ < P2 < Ps, we can show that 
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-1 < yU < 0. This constraint on ju is all that is needed to derive the main results of this paper, so the precise placement 
of the surface layer is not important for our argument below. On the other hand, moving the potential vorticity layers 
around would necessitate non-symmetric generalizations of the operator ^ap, which may be interesting for oceano- 
graphic modeling, but not necessary for atmospheric modeling, and beyond the scope of this paper. We will therefore 
assume that p\ and p2 are fixed but allow to vary as pi < Ps < lAtm, which in turn corresponds to -1/3 < ;U < 0. 

The dissipation term configuration given by Eq. (ISFt and Eq. ( |59] l corresponds to setting the generalized dissipation 
operator spectrum Dap{k) equal to 



D(k) 



Di(k) 
lud(k) Diik) + d(k) 



with Di(q), D2(q), and d(q) given by 

Di(k) = vk^P^^ and D2(k) = (v + Av)k^P^^ and d(k) = Avsk^. 
Note that for jj -0 and A - 1, this reduces to the simpler case of streamfunction dissipation. 



(61) 



(62) 



4.2. Dissipation rate spectra for the two-layer model 

We may now leverage Eq. ( |39] l and Eq. ( l40l l to calculate the energy and potential enstrophy dissipation rate 
spectrum DE{k) and Dc{k) in terms of the streamfunction spectra Ui{k), Utik), and Cnik). For the case of the energy 
dissipation rate spectrum Dgik), noting that Dnik) = 0, a simple calculation gives 



DE{k) = 2Dn(k)Ui(k) + 2D22(k)U2(k) + 2D2i(k)C2i(k) 

= 2Dn(k)Ui(k) + 2D22(k)U2(k) + D2i(k)[2C nik) - U(k)] + D2i(k)U(k) 

= [2Dn(k) + D2x{k)]U,{k) + [2D22(k) + D2i(k)]U2(k) + D2i(k)[2C n(k) - U(k)] 



= Af(k)Udk)+A'^'(k)U2(k)+A'^'(k)[2Ci2ik) - U(k)] 

with A^^\k), Af{k\ and Af(k) given by 

A^^\k) = 2Dn{k) + D2i{k) = 2Diik) + ndik), 
Af(k) = 2D22{k) + D2i(k) = 2D2{k) + 2d{k) + pd{k), 
Af{k) ^ D2i(k) ^ Ijd(k). 



l(2) 



(63) 
(64) 
(65) 
(66) 



(67) 
(68) 
(69) 



We note that terms involving the streamfunction cross-spectrum Cnik) have been reorganized in terms of 2Ci2ik) - 
U{k) so that we can take advantage of the inequality 2Ci2(fc) - U{k) < 0. For the potential enstrophy dissipation rate 
spectrum D(j{k), we take advantage of the symmetry assumption Lapik) - L^aik) to rewrite Eq. (l40l l as 



Dcik) = -2Y^L,p{k)D„y{k)Cpy{k) = -2Y^Lp,{k)Day{k)Cpy{k) = -2Yj.LD)py{k)Cpy{k). 



The components of (LD)(k) are given by 



(LD)(k) 



a(k) 
b(k) 



b(k) 
a(k) 



Di(k) 
jidik) 



a(k)Di(k) + nb(k)d(k) 
b{k)Dx(k) + na{k)d{k) 



and it follows that Dc{k) is given by 





D2(k) + dik) 

b(k)[D2(k) + d(k)] 
a(k)[D2(k) + d(k)] 



Dcik) = -2{iLD)nik)Uiik) + (LD)22(k)U2(k) + [(LD)n{k) + (LD)2i(k)]Cn(k)] 

= -{[2(LD)n(k) + iLD)n(k) + (LD)2i(k)]Ui(k) + [2(LD)22(k) + {LD)n{k) + {LD)2i(k)]U2{k) 
+ {{LD)n{k) + {LD)2i{kmCi2{k) - U{k)]] 



(70) 

(71) 
(72) 

(73) 
(74) 
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= A-^\k)Um + A'^'ikWiik) + A'^\k)[2Cn{k) - U{k)l (75) 
with A^^\k), A^^\k), A^^\k) given by 

Al\k) = -Yl(LD)u(k) + (LD)n(k) + (LD)2x{k)] (76) 

= 2[a{k)D]ik) + i^b(k)d(k)] + b(k)Di(k) + i^a(k)d(k) + b(k)[D2(k) + d(k)] (77) 

= [2a(fe) + b(k)]Di (k) + b(k)D2(k) + [2b(k) + a{k)]ixd{k) + b{k)d{k\ (78) 

Ag)(A:) = -[2{LD)22{k) + {LD)n{k) + {LD)2i{k)] (79) 

= 2a{k)[D2{k) + d{k)\ + b{k)Dy{k) + na{k)d{k) + b{k)[D2{k) + d{k)\ (80) 

= bik)Diik) + [2aik) + b{k)\D2{k) + [2a{k) + b{k)\d{k) + na{k)d{k), (81) 

Afik) = -[iLD)uik) + (LD)2i(A:)] = b{k)Di{k) + iia{k)d{k) + b{k)[D2{k) + d{k)\ (82) 

= b{k)[Di{k) + D2ik)] + b{k)d{k) + iia{k)d{k). (83) 



The above expressions for Dsik) and Dc(k) are the point of departure for the investigation of the flux inequality 
under the general case of streamfunction dissipation with extrapolated Ekman dissipation and difi'erential small-scale 
dissipation. 

4.3. Sufficient conditions in terms of dissipation coefficients 

As we have discussed previously, to satisfy the flux inequality k^IlEik) - Tldk) < for a given wavenumber k, 
it is sufficient to show that A{k, q) < for all wavenumbers k < q. Using our previous expressions for the energy 
dissipation rate Dsik) and the potential enstrophy dissipation rate Ddk), we can calculate A{k, q). Consequently, 



A{k, q) is given by 

Aik,q) = k^DE{q)-DG(q) (84) 

= Ai(k, q)Ui(q) + A2(k, q)U2(q) + Ai(k, q)[2Cn(q) - U(q)], (85) 

with Ai{k, q), A2{k, q), and AT,{k, q) given by 

A,{k,q) = eA^^\q)-Al\q) (86) 

= e\_2D, (q) + fxdiq)] - [2aiq) + b{q)]DM - b{q)D2{q) - \_2b{q) + a{q)]nd{q) - b{q)d{q) (87) 

= [2^ - 2a(q) - b{q)-\DM) - Kq)D2{q) - b{q)d{q) + - 2b{q) - a{q)\nd{q\ (88) 

A2{k,q)^l^Af{q)-A'^\q) (89) 
= k\2D2{q) + Hd(q) + 2d(q)] - b(q)Di(q) - [2a(q) + b(q)]D2(q) 

- [2a(q) + b(q)]d(q) - ^ia(q)d{q) (90) 

= -b(q)Di(q) + \2k^ - 2a(q) - b{q)]D2{q) + \2k^ - 2a(q) - b(q)]d(q) + ii[k^ - a{q)]d{q), (91) 

A^{k, q) = k^Afiq) - A'^\q) = k^iidiq) - b(q)[Di(q) + D2(q)] - b{q)d(q) - na{q)d(q) (92) 

= -biq)[Diiq) + D2m - b{q)d{q) + fiik" - a{q)]d{q). (93) 



We observe that Ui{q) > and U2{q) > and 2Cniq) - U{q) < 0, consequently the sign of A{k, q) depends on the 
sign of the coefficients A\{k, q), A2ik, q), and A^ik, q). For the argument below we may assume that -1 < // < and 
D\(q) < D2(q). Here Di(q) < D2(q) corresponds to differential small-scale diffusion (i.e Av > 0) and Di(q) = D2(q) 
corresponds to symmetric small-scale dissipation (i.e Av - 0). We begin our argument with the following lemma, 
stating essentially that A2ik, q) is negative, as long as Di{q) < D2{q), and A3 (A:, q) is always positive. 

Lemma 1. Assume that b{q) < and k^ - a{q) - b{q) < 0. Assume also streamfunction dissipation with both 
differential small-scale dissipation and extrapolated Ekman dissipation with —1 < p < Q. Then A-i{k,q) > 0, and 
furthermore, ifDi(q) < D2(q), then we also have A2(k, q) < 0. 
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Proof. We recall that A^ik, q) is given by 



A^ik, q) = -b{q){Di{q) + D2{q)] - b(q)d(q) + n{k^ - a{q)]d{q). (94) 

Since, by definition, D\(q) > 0, and D2{q) > 0, and d{q) > 0, and since b{q) < 0, and k^ - a(q) - [k^ - a(q) - b(q)] + 
b(q) < k^ - a(q) - b(q) < 0, and fi < 0, it follows that all contributing terms to A^ik, q) are positive and therefore 
AT,{k, q) > 0. For the case of A2{k, q), let us assume first that D\{q) < Diiq). We rewrite A2{k, q) as follows: 

A2{k, q) = -b(q)Di(q) + [2k^ - 2a(q) - b{q)]D2(q) + [21^ - 2a{q) - b(q)]d(q) + fi[k^ - a(q)]d(q) (95) 
= -b(q)[Di(q) - D2(q)] + 2[k^ - a(q) - b(q)]D2iq) 

+ 2[k^ - a(q) - b(q)]d(q) + b(q)d(q) + ^i[k^ - a(q) - b(q)]d(q) + ^ib(q)d(q) (96) 

= -b(q)[D,(q) - D2{q)] + - a(q) - b(q)][2D2{q) + (2 + ^i)dm + (M + mq)d(q). (97) 

Since b{q) < 0, and k^ - a{q) - b(q) < 0, and yu + 2 > 0, and fi + I > 0, we see that all contributing terms to A2{k, q) 
are negative and therefore A2{k, q) < 0. This concludes the proof. □ 

Proposition 2. Assume streamfunction dissipation with both differential small-scale dissipation and extrapolated 
Ekman dissipation with -1 < < 0. Assume also that k^ — a{q) — b{q) < 0, and b{q) < 0, and Diiq) > 0, and 
D2{q) > 0, and AD{q) = D2{q) - D\{q) > 0, and also that D\{q), AD(q), and d(q) satisfy 

2D,(q) + fld(q) ^ biq) 
AD(q) + + \)d{q) k^ - a(q) - b(q) ' 

Then it follows that A(k, q) < 0. 

Proof. We recall that A(k, q) is given by 

Aik, q) =Adk, qWiiq) + A2(k, q)U2(q) + Aiik, q)[2Cn(q) - Uiq)]. (99) 

Using the previous lemma, from the given assumptions above, we have A2(k,q) < and Ai(k,q) > 0. Now let us 
rewrite A i {k, q) as 

Ai(k, q) = [2k^ - 2a(q) - b(q)]Di{q) - b(q)D2(q) - b{q)d(q) + [k^ - 2b(q) - a(q)]iid{q) (100) 
= 2[k^ - a(q) - b(q)]Di(q) - b(q)[D2(q) - Di{q)] - b{q)d{q) 

+ [k^ -a(q)- b{q)]nd{q) - iub(q)d(q) (101) 

^ [^2 _ a(q) - b(q)][2Di(q) + jidiq)] - b(q)[AD(q) + (fi + l)d(q)]. (102) 

Since k^ - a(q) - b(q) < and AD(q) + (fi + l)d(q) > 0, it follows that A i (A:, ^) < if and only if Eq. (|98] l is satisfied. 
Thus, since we also know that Ui(q) > 0, and U2(q) > 0, and 2Ci2(^) - U{q) < 0, it follows that all terms contributing 
to A(k, q) are negative, and therefore A{k, q) <Q. □ 

It is easy to see that using Di(q) - vq^^^^, and d{q) - {3/2)vEq^ (the 3/2 factor follows from using A - 3/2), and 
AD{q) = D2{q) - Di{q) - Avq^P^^, and fi - -1/3, proposition|2]gives the following statement 

Avq^P + Vf q^ - k^ 

^<t\ -<^^—^^A{k,q)<Q. (103) 

Avq^P - ve kj^ 

Note that the hypothesis given by Eq. ( I103l l requires that ve < Avq^'', which allows us to invert both sides of Eq. ( I103l l 
to yield Eq. ( |98] ), so the proposition can be then applied. For ve > 4-vq^P, the Eq. ( |98] l is violated, as the left-hand 
side becomes negative while the right-hand side remains positive. More precisely, in Eq. (|98T l, the right-hand side is 
positive for q > k, the denominator of the left-hand side satisfies AD{q) + (fJ. + l)d(q) > by the given hypotheses, 
and the constraint ve < A-vq'^P is needed to ensure that the numerator 2D\{q) + ^d{q) is not negative, so that it can 
be possible for Eq. ( |98] | to be satisfied. Consequently, we see that increasing either ve or Av indicates a tendency 
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towards violating the flux inequality. The role of differential diffusion is very important here since, for Av > 0, the 
left-hand-side of the hypothesis in Eq. (1103b will approach Av/(4v) and remain bounded for large wavenumbers q, 
whereas for Av - 0, the same left-hand-side will vanish rapidly with increasing wavenumber q. As a result, violating 
the flux inequality may become easier under differential small-scale dissipation Av. On the other hand, the role of 
ve becomes even more dramatic, since increasing ve from towards Avq^'' will result in a hyperbolic blow-up of the 



left-hand-side of the hypothesis of Eq. ( 1103b . thus yielding an even more rapid violation of the hypothesis. 

The presence of ve in the denominator of the left-hand-side fraction of Eq. (1103b is due to the use of extrapolated 
Ekman dissipation. For the case ^-Q and A-\of streamfunction dissipation, proposition |2]reduces to the statement 
given by 

Avq^P + vf q^ - 

V^'^''^^'"' ''''' 

where the hyperbolic blow-up is no longer possible. Differential small-scale dissipation however maintains its ten- 



dency towards violating the flux inequality for increasing Av since the left-hand-side in the hypothesis of Eq. ( 1104b 
still approaches Av/(4v) in the limit of large wavenumbers q. 

It is also interesting to consider yet another special case in which we eliminate differential small-scale dissipation 
but retain extrapolated Ekman dissipation. This corresponds to choosing fi = -1/3 and Av = 0, and the statement of 
Eq. ( 11031 ) can be now simplified to read 

''''' 

Now let us compare that against the case where both differential small-scale dissipation and extrapolated Ekman 
dissipation are eliminated (i.e. fi - and Av = 0). The corresponding sufficient condition is given by 

Ve q^ - 

We see that the right-hand-side in the hypothesis of Eq. (11061 ) grows quadratically with increasing wavenumber q, 
thereby making it very easy to satisfy the hypothesis regardless of how much one increases the Ekman dissipation co- 
efficient Ve- On the other hand, under extrapolated Ekman dissipation, even without help from differential small-scale 



dissipation, the right-hand-side of the hypothesis of Eq. ( 1105b remains bounded in the limit of increasing wavenum- 



ber q. Consequently, it is clear that for sufficiently large V£, the hypothesis of Eq. ( 1105b can be easily violated, and 
including differential small-scale dissipation should make it all the more easier 

It should be stressed that in the above discussion, the hypotheses given by Eq. ( 1103b . Eq. ( 1104b . and Eq. ( 1105b 
are sufficient conditions, however they are not necessary conditions. A violation of Eq. ( |98] | will ensure that the 
term A\{k,q)U\(q) gives a positive contribution to A{k,q). However, according to Lemma [T] the contributions of 
A2{k, q)U2{q) and A3(A;, q)[2Ci2{q)-U{q)] will remain negative, so the sign of A{k, q) is dependant on which term gets 
to be dominant. Therefore, it is far from a foregone conclusion that a violation of the flux inequality is possible under 
the dissipation configurations considered above. However, the significant tightening of the sufficient condition with 
the introduction of extrapolated Ekman dissipation and differential small-scale dissipation indicates that a violation of 
the flux inequality may be becoming easier to achieve, under these configurations. 



4.4. Sufficient conditions in terms of streamfunction spectra 

We would now like to consider statements providing sufficient conditions for satisfying the flux inequality, for- 
mulated in terms of the streamfunction spectra U\(q), U2(q), and Cniq), for the dissipation configuration given by 
Eq. (|58] | and Eq. (|59] l, i.e. streamfunction dissipation with differential small-scale dissipation and extrapolated Ekman 
dissipation. Such conditions imply corresponding constraints on the distribution of energy and potential enstrophy 
between layers, to be explored in future work. We will derive propositions for three separate cases. Proposition |3]cor- 
responds to streamfunction dissipation without differential small-scale dissipation and without extrapolated Ekman 
dissipation. Proposition|4]corresponds to streamfunction dissipation with both differential small-scale dissipation and 



15 



extrapolated Ekman dissipation. Finally, Proposition |5] corresponds to streamfunction dissipation with extrapolated 
Ekman dissipation but without differential small-scale dissipation. 

The first step towards deriving the propositions below is to rewrite A{k,q) in terms of Di(q), D2{q), and d{q) as 
follows 

A(k,q) = Bi(k,q)Di(q) + B2(k,q)D2(q) + B^{k,q)d{q), (107) 

with B\{k, q), B2{k, q), and B^ik, q) given by 

Bi(k,q) = [2^ - 2a(q) - b(q)]Ui(q) - b(q)U2(q) - b(q)[2Cn{q) " U{q)] (108) 

= 2[k^ - a(q)]Uiiq) - 2b{q)Cn{q\ (109) 

B2(k, q) = -b(q)Ui(q) + [2^ - 2a(q) - b(q)]U2(q) - b(q)[2Ci2(q) - U(q)] (110) 

= 2[k^ - a(q)]U2(q) - 2b(q)Cn(q), (HI) 
B^ik, q) = -biq)Ui(q) + [k^ " a(q) - 2b(q)]^iU i(q) + [2k^ - 2a{q) - b(q)]U2(q) 

+ flik" - a(q)]U2(q) - b(q)[2Cn(.q) ' U(q)] + /.[^ - a(qmC niq) ' U(q)] (112) 
^2[k^-a(q)]U2(q)-2b(q)Cn(q) 

+ fi[k^ - a(q)][Uiiq) + Uiiq) + 2Cn(q) - U(q)] - 2fibiq)U,iq) (113) 

= 2[k^ - a(q)]U2(q) - 2biq)Cn(q) + tiik" - a{q)]2Cn{q) - 2^ibiq)U liq) . (114) 

More specifically, we were able to derive the following three propositions from the above equations. 

Propositions. Assume symmetric small-scale streamfunction dissipation with d(k) > and D\{k) — D2(k) = D(k) > 
0, standard Ekman term (i.e. - 0) and k^ - a{q) — b(q) < and b{q) < 0. It follows that ifCniq) ^ U2{q), then 
A(k,q) < 0. 

Proof We write A(k, q), under the assumption of symmetric small-scale dissipation (i.e. Di (q) = D2iq) = D(k)), as 

Mk, q) = [BiiK q) + B2(k, q)]D(q) + B^ik, q)d{q). (115) 

We note that from the given assumptions, we have 

Bi{k,q) + B2(k,q) = 2[k^ - a{qWi{q) - 2b(q)Cniq) + 2[^' - a(q)]U2(q) - 2biq)Cn(q) (116) 
= 2[k^ - a(q)]U(q) - 4b(q)Cn(q) (117) 
= 2[k^ - a(q) - b{q)]U{q) - 2b{q)[2Cn{q) - U{q)] < 0, (118) 

using k^ - aiq) - b(q) < 0, b(q) < 0, and 2C\2{q) - U(q) < 0. From the hypothesis Ci2iq) ^ Ui^q), we can also show 
that 

B^iKq) = 2[e - a{q)]U2{q) - 2b(q)Cn(q) (119) 
< 2[k^ - a(q)]U2(q) - 2b{q)U2(q) = 2[k^ - a(q) - b{q)]U2{q) < 0. (120) 

Since d{q) > and D{q) > 0, it follows that A{k, q) <0. □ 

This proposition shows that under symmetric small-scale streamfunction dissipation alone, using standard as op- 
posed to extrapolated Ekman dissipation, the inequality Cniq) < U2{q) implies A(k, q) < for all wavenumbers 
k < q. We already know that Cniq) is mathematically restricted via the triangle inequality 2\Ci2(q)\ < Ui(q) + U2{q) 
over an interval of values intersecting with the constraint Cniq) < U2iq), so the actual constraint on Ci2(<?) is tighter. 

We were able to find similar sufficient conditions for the cases of extrapolated Ekman dissipation and differential 
small-scale dissipation, however including either or both features on top of streamfunction dissipation makes the 
sufficient conditions more restrictive. This is, of course, expected and consistent with the preceding discussion of the 
consequences of Proposition |2l The corresponding propositions and their proofs read: 
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Proposition 4. Assume that b{q) < and - a{q) — b{q) < 0. Assume also streamfunction dissipation with both 
differential small-scale dissipation and extrapolated Ekman dissipation with -1 < yu < 0. 

1. IfCniq) < 0, then A(k,q) < 0. 

2. IfCniq) < min{t/i(^), U2(q)] and Uiiq) + jiUiiq) > 0, then Aik, q) < 0. 

Proof. To show (a) we first note that - a(q) = [k^ - a(q) - b(q)] + b{q) <k^ — a(q) - b{q) < 0. Combined with the 



given assumptions, we find that Bi(k, q), B2(k, q) and BT,{k, q) satisfy 

B^{k,q) - 2[k^ - a{qWM) ' 2b{q)Cn(q) < 2[k^ - a(q)]Ui(q) < 0, (121) 

B2(k,q) = 2[fc2 - a(q)]U2iq) " 2b(q)Ct2iq) <2[k^- a(q)]U2(q) < 0, (122) 

B3ik,q) = 2[k^ - aiqWiiq) " 2b{q)Cx2{q) + l^[k^ - a{q)\2Cniq) - 2fibiq)Uiiq) (123) 

< 2[k^ - a(q)Wi{q) - 2iib{q)Ui{q) < 0. (124) 

Here we used the inequaUties -2b(q)Ci2(q) < 0, and iJ.[k^ - a{q)]2Ci2iq) < 0, and 2ij.b{q)Ui{q) < 0, that follow from 
the given assumptions. It follows that A(k, q) < 0. 

To show (b) we use the given assumptions to show that 

Bi(k, q) = 2[k' - a{q)WM - 2biq)Cuiq) < 2[fc^ - aiq)]Uiiq) - 2b{q)Ui{q) (125) 

< 2[1^ - a{q) - b{q)WM) ^ 0, (126) 
B2{k, q) = 2[k' - aiq)]U2iq) - 2biq)Cuiq) < 2[fe^ - aiq)]U2iq) - 2b{q)U2{q) (127) 

= 2[A^ - a{q) - b{q)W2{q) < 0. (128) 

The above two inequalities for B\{k, q) and B2{k, q) are based on the assumptions Cniq) < U\{q) and Cniq) ^ Uiiq)- 
We also show that B3{k, q) is bounded by 

B^iKq) = 2[k' - aiq)]U2iq) - 2b{q)Cn{q) + tAk" - aiq)]2Cniq) - 2ixbiq)Uiiq) (129) 

< 2[k' - aiq)]U2iq) - 2b{q)U2{q) + fxik" - a{q)\2Ui{q) - 2nbiq)Uiiq) (130) 
= 2[k^ - a{q) - biq)]U2iq) + 2iAk^ - a{q) - b{q)Wi{q) (131) 
= 2\e - a{q) - b(q)][U2(q)+fiUim- (132) 

Since k^ - a{q) - b{q) < and by hypothesis Uiiq) + ixU2iq) > it follows that B3{k,q) < 0, and consequently 

A{k,q)<Q. □ 



Proposition 5. Assume that k^ — a(q) — b(q) < and b(q) < 0. We also assume streamfunction dissipation with 
extrapolated Ekman dissipation with -1 < /i < and symmetric small-scale dissipation with Di{q) - D2{q) > 0. It 
follows that ifCniq) < min{L''i(^), U2iq)} then A{k, q) < 0. 

Proof Under the assumption of symmetric small-scale dissipation (i.e. D\{q) = D2{q)), we may rewrite A{k, q) as 



A{k, q) = [By{k, q) + B2(k, q)]D(q) + Bjik, q)d{q). (133) 

We note that from the given assumptions, we have 

Bi{k,q) + B2{k,q) = 2[k' - a(q)]Ui{q) - 2b{q)Cn(q) + 2[k' - a{q)W2{q) - 2b{q)Cn{q) (134) 
= 2[k^ - a(q)]U{q) - 4M<?)Ci2(?) (135) 
= 2{k^ - a{q) - b{q)]U{q) - 2b{q)[2Cn{q) - U{q)] < 0, (136) 

using A:^ - a{q) - b{q) < 0, b{q) < 0, and 2Cu{q) -U{q) < 0. We also have 

B^ik, q) = 2[k' - aiq)]U2iq) - 2biq)Cuiq) + pik" - aiq)]2Cuiq) - 2pb{q)UM) (137) 
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< 2{k^ - aiqW2iq) - 2b(q)U2(q) + ^i[k^ - a(q)]2Ci2(q) - 2nb(q)Cn(q) (138) 
= 2[k^ - aiq) - b(q)]U2(q) + 2fi[k^ - a(q) - b(q)]Cn(q) (139) 

< 2[k^ - a(q) - b(q)]U2(q) + 2f,[k^ - a(q) - b(q)]U2(q) (140) 
= 2(1 + ,,)[k^ - a(q) - b(q)]U2(q) < 0. (141) 

Here, on the first line we used the assumptions Cniq) < U\{q) and Cniq) < U2(q) to argue that -2biq)C\2{q) < 
-2b(q)U2{q) and -2iJ.b(q)U\(q) < -2fib(q)Ci2{q)- The remainder of the argument continues to apply the given 
assumptions and it is easy to follow. We conclude that A{k, q) <Q. □ 



Proposition|4]shows that for the more general case of streamfunction dissipation with both differential small-scale 
dissipation and extrapolated Ekman dissipation, if the streamfunction spectrum Cniq) is negative for all wavenumbers 
q > k, then the flux inequality is satisfied at the wavenumber k. It also shows that the restriction on the streamfunction 
spectrum Cniq) can be relaxed to the wider inequality Cniq) ^ min{f/i(^), U2iq)} if we choose to introduce the re- 
striction Uiiq)+ fiU2iq) > on the streamfunction spectra Uiiq) and U2iq)- In proposition|5]we eliminate differential 
small-scale dissipation but retain extrapolated Ekman dissipation. This eliminates the restriction Uiiq) + fiU2iq) ^ 
whereas the restriction on the streamfunction spectrum Cniq) remains the same as in Proposition]?] This shows that 
the restriction U\iq) + fiU2iq) > originates from differential small-scale dissipation, and since fi is negative, it con- 
stitutes a non-trivial constraint on the streamfunction spectra Uiiq) and U2iq)- As a result, the sufficient conditions 
of Proposition|5]are rigorously wider than the sufficient conditions of Proposition]?] It goes without saying that elim- 
inating both differential small-scale dissipation and extrapolated Ekman dissipation reverts us back to Proposition |3] 
where the stated sufficient condition is clearly wider than that of Proposition |5] Specifically, for /v = 0, the inequality 
U\iq) + iiU2iq) > reduces to the trivial inequality U\iq) > 0. Furthermore, in the proof of proposition |5] if ju = 0, 
we no longer need the constraint Cniq) ^ Uiiq) to show that B^ik, q) < 0, and only the constraint Cniq) ^ U2iq) is 
needed by the remainder of the proof. 

5. Conclusions and Discussion 

We have derived rigorous sufficient conditions for satisfying the flux inequality k^YlEik) - Tldk) < for a general 
«-layer quasi-geostrophic model with constant layer-by-layer thickness, under symmetric streamfunction dissipation. 
By symmetric streamfunction dissipation we mean that for every layer the dissipation term is given by the same linear 
Fourier-diagonal operator, applied only on the streamfunction field of the same layer. It follows that under symmetric 
configurations of the dissipation terms, n-layer quasi-geostrophic models will indeed have a phenomenology similar 
to two-dimensional Navier-Stokes turbulence. Asymmetric dissipation configurations, where different dissipation op- 
erators are used on different layers, have been considered for the special case of a two-layer quasi-geostrophic model, 
dissipated with general streamfunction dissipation with or without extrapolated Ekman dissipation and differential 
small-scale dissipation. We have derived results of the form that if the degree of asymmetry in the dissipation terms 
between the two layers is bounded to a certain extent, then the flux inequality will continue to be satisfied. Our results 
on the non-trivial dependence of the dissipation rate spectra of energy and potential enstrophy on the energy and 
potential enstrophy spectra via the streamfunction spectra, are also very relevant to the correct formulation of closure 
models for multi-layer quasi-geostrophic systems. 

One limitation of the current investigation is that we have disregarded the beta term, mainly to avoid the math- 
ematical difficulties associated with the anisotropic nature of the term. This elimination can be tolerated, from a 
physical standpoint, as long as the beta term is active only in the forcing range and the baroclinic forcing at the same 
forcing range is powerful enough to overshadow the beta term. As long as the effect of the beta term remains limited 
to large scales (i.e. planetary and synoptic scales), it will not contribute to the integrals of Eq. (l35T l and Eq. (l36l l and 
the results reported in this paper will remain entirely unaffected. The only other physical assumption inherent in these 
results is that forcing via the baroclinic instability is limited to large scales. The propositions 1-5 are mathematically 
rigorous and do not require these assumptions, however the assumptions come into play at the very last step where 
the conclusion of propositions 1-5 (i.e. Aik,q) < 0) is used to infer the flux inequality itself. On the other hand, no 
phenomenological assumptions about any spectrum are needed at any step of the argument. 
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For the case of the two-layer quasi-geostrophic model, we have seen that, starting from a streamfunction dissi- 
pation configuration, adding either extrapolated Ekman dissipation or differential small-scale dissipation (or both) 
tends to tighten the sufficient conditions for satisfying the flux inequality. This indicates that the flux inequality 
k^TlEik) - Ilcik) < may be more easily violated under these more general dissipation configurations. A violation 
of the flux inequality beyond a wavenumber k, would then allow a downscale energy flux larg e enough to result in a 
transition from k^^ to k^^^^ scaling in the energy spectrum near the wavenumber k, Lioi II, 3^. However, while there 



is a plausible physical motivation for using extrapolated Ekman dissipation, there is no obvious physical motivation 
for introducing an asymmetric configuration of the small-scale dissipation terms. We would therefore like to expand 
on the reasons why we believe that this is an idea worth pursuing. 

Any kind of small-scale dissipation in quasi-geostrophic models is not physical but is tolerated mainly because it 
is intended to model the dissipative mechanisms that exist at smaller scales where quasi-geostrophic dynamics breaks 
down and three-dimensional dynamics becomes dominant. Lindborg f3l\ estimates that quasi-geostrophic dynamics 
break down at a length scale of about 100km. However, the scaling transition wavenumber kt of the Nastrom-Gage 
spectrum [16-19], and consequently the breakdown of the flux inequality, occurs at a greater length scale of about 
1000km to 700km in wavelength, which is still within the quasi-geostrophic regime. The hypothesis underlying the 
quasi-geostrophic modeling of atmospheric turbulence is that the locality of the coexisting downscale potential en- 
strophy cascade and downscale energy cascade shields them from the three-dimensional dominated regime at the 
smallest scales. Both cascades are furthermore protected by the continuing conservation of potential enstrophy under 
the stratified turbulence dynamics that becomes dominant at scales less than 100km. The above considerations suggest 
the hypothesis that three-dimensional effects will not contaminate the nonlinear quasi-geostrophic dynamics driving 
the coexisting cascades of potential enstrophy and energy in the quasi-geostrophic regime, which allows us to model 
smaU-scale three-dimensional processes, as seen from the quasi-geostrophic regime's point of view, via small-scale 
hyperdiffusion terms applied to all layers. That said, there is the non-local effect that the anomalous energy sink, pro- 
vided by the three-dimensional regime, can inflict on the quasi-geostrophic regime, and that is to boost the downscale 
energy dissipation rate, thereby increasing the downscale energy flux passing through the quasi-geostrophic regime 
of length scales and moving the transition wavenumber kt deep into the inertial range. This is why we propose that if 
differential small-scale dissipation can be shown to achieve an equivalent effect, then it should be accepted as a more 
realistic configuration, for modelling purposes. 
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Appendix A. Proof that E(k) is always positive 

In this appendix we prove the inequality 2\Cafi{k)\ < Ua(k) + Up{k), and we use it to show that if the matrix L^pik) 
satisfies the following conditions: 

> 0, for a /3, (A.l) 
Yj^ap{k)<Q, (A.2) 

then the energy spectrum E{k) will be always positive with E{k) > 0. 
To establish the inequality, we first note that 

(lA,, ± i/r,, + lA/j)^ = <iA«, ^a)k ± 2 (lA^, + ("A/J, = Ua(k) + U^ik) + ICapik). (A.3) 

Since {ij/a ± if/p, + ij/p)^ > 0, it follows that Ua(k) + Upik) - 2Cap{k) > and Ua(k) + Upik) + 2Cap(k) > 0, and 
therefore we have 



2\Cap{k)\ < Ua(k) + Up{k). 
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(A.4) 



Now, to show that the energy spectrum E(k) is always positive, we begin by rewriting Eq. (129) as follows: 

Eik) = - ^ Up{k)Cap{k) = - 2 Ua(k)Uaik) - 2 Up{k)Cap{k) (A.5) 



= - ^ LUk)Ua(k) - i ^ L,p(k)[Ua(k) + Up(k)] - i ^ L,^(^)[2C„/j(^) - UAk) - Up(k)] (A.6) 
= - ^ LUkWaik) - i + Lp,(k)]U,ik) - i ^ L„/5(^)[2C„/5(^) - U^k) - Upik)] (A.7) 

ap ap 

= - ^k»(^) + ^ + Lp„{k)-h,{k) - i 2 L„^(fe)[2C„^(^) - - Upik)] (A.8) 

tr /3 ap 

a+P a+p 

= - [Z ^"'^(^^l ^"^"^^ -\Yj ^-pWCapik) - Ua(k) - Up{k)] (A.9) 

a p ap 

We note that Ua{k) > 0, since Ua(k) is always positive, and we have just shown that 2Cap(k) - Ua(k) - Upik) < 0. 
Combining these with the assumptions given by Eq. ( lA.lb and Eq. ( IA.2b . we see that both terms in our expression for 
E{k) are positive and therefore E{k) > 0. 

Appendix B. Derivation of dissipation rate spectra 

In this appendix, we will show that the energy dissipation rate spectrum Dsik) and the layer-by-layer potential 
enstrophy dissipation rate spectra Dc^^ (k) are given by 

DE(k)^2Y,D„p{k)Cap{k), (B.l) 

ap 

Dc„ (k) = -2 Z L„p(k)Da,y{k)Cpy(k). (B.2) 

Py 

The proof mirrors the argument used in my previous paper [5o^ to derive the energy forcing spectrum and the potential 
enstrophy forcing spectrum for the same model. We begin by writing the governing equation for the streamfunction 
field 41 a as 

^ + Z ^apJi^P^ qp) = ^^p ^py'^y + Z ^»pfp- ^^-^^ 

P Py P 

Differentiating the streamfunction spectrum Capik) with respect to time gives 
dCgpjk) _ I dip, \ I dM 



8t \5.'H"r'^/.' ^^-'^ 

and we may write a governing equation for Capik) in the form: 
dCapik) 



dt 



+ 7apik) = -'Dapik) + y^pik). (B.5) 



Here, 7apik) is the contribution from the nonlinear Jacobian term, D^pik) is the contribution from the dissipation term, 
and "Japik) is the contribution from the forcing term. The dissipation term Dapik) can now be obtained by replacing 
in Eq. ( IB.4I ) the streamfunction time-derivative dij/aldt with the dissipation term ^pyifiy. This gives 

•Dapik) = - (Z -^^y ^ys^s, M - Ua, J] ^Py ^r^Ai ) (B.6) 
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= - Y}l^a\ik)Dys{k)Cpm + L-^l(,k)Dys{k)Ca5m- (B.7) 
js 

We may now easily write the dissipation rate spectra Dsik) and Dc{k) by applying on Dapik) the linear operators 
indicated by Eq. ( l29b and Eq. ( l32b . We therefore find that the energy dissipation rate energy spectrum Dsik) is given 
by 

DE(k) = - 2 La4k)V„fiik) = 2 [L„^(fe)L-;(fe)D^,(^)C;,5(fe) + L„/j(fe)L^;(fe)D^i(^)C„a(fe)] (B.8) 

apy6 

" Z[Z ^'/^«(^)'L„;(^)]£>r*(^)C/s^(^) + X[Z i'.-/j(^)'L^r(^)]'Dri(^)C,<5(^) (B.9) 

= J] (5;j^Dy5(A:)C/j5(A:) + J] 5,^Z)^5(A:)C<,i(fc) (B. 10) 

= ^ Dps(k)Cps{k) + 2 D^i(^)Cya(^) = 2 ^ D„;s(fc)C„^(^). (B.l 1) 

/SS yd a/S 

The layer-by-layer potential enstrophy spectrum Dg„ (fe) is likewise given by 

DcAk) = 2 L,p{k)Uy{k)'Dpy{k) = - 2 Up{k)Uy{k)[L-pl(k)Dse{k)Cy,{k) + L;](A:)D&(A:)C;,,(^)] (B.12) 

= - 2 <J.rfi«r(fc)£>a.(fc)C^.(^) - J] 5„5V(^)Z)&(^)C;j,(fe) (B.13) 

ydE /3S£ 

= - 2 L„y(k)DUk)CyM - ^ L„is(k)DUk)Cf}M (B. 14) 

= -2 ^ Up{k)D„y{k)Cpy{k). (B. 15) 

The corresponding conservation laws read 
dE{k) dllEik) 



dt dt 
dG{k) dlicik) 

dt dt 



DE(k) + FE(k), (B.16) 
Dc(k) + Fc(k). (B.17) 



We see that positive Dsik) and Dc(k) correspond to the case where the dissipation terms are truly dissipative. This 
concludes the argument. 
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